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Abstract 

We present a first calculation of the generalized parton distributions of the photon (both polarized 
and unpolarized) using overlaps of light-front wave functions at leading order in a and zeroth order 
in as', for non-zero transverse momentum transfer and zero skewness. We present the novel parton 
content of the photon in transverse position space. 
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Introduction 



In deep inelastic scattering of a highly virtual photon on a real photon, the partonic con- 
stituents of the photon play a dominant role when the virtuality is very large. In this 
case the pointlike contribution to the photon structure function F2{x,Q'^) dominates over the 
hadronic component. This can be calculated perturbatively. These are relevant in the context 
of e^e~ annihilation and photoproduction. Unlike the proton structure function, where only 
the dependence is calculated perturbatively and the x dependence has to be fitted using 
experimental data, in the photon structure function, both x and dependence can be calcu- 
lated. F2{x,Q'^) show logarithmic dependence already in parton model, unlike the proton 
structure function Leading order QCD calculation differs from the parton model result for 
F2{x,Q'^) by calculable finite terms |2|. The photon structure function is now known fairly 



accurately and agrees well with experimental results 

In j4| deeply virtual Compton scattering (DVCS) 7*7 — )■ 77 on a photon target was con- 
sidered in the kinematic region of large center-of-mass energy, large virtuality (Q^) but small 
squared momentum transfer {—t). The result was interpreted at leading logarithmic order as a 
factorized form of the scattering amplitude in terms of a hard handbag diagram and the gener- 
alized parton distributions of the photon. The calculation was done at leading order in a and 
zeroth order in as when the momentum transfer was purely in the longitudinal direction. These 
are called anomalous GPDs as they show logarithmic scale dependence already in parton model. 
They are of particular interest as they can be calculated in perturbation theory and can act as 
theoretical tools to understand the basic properties of GPDs like polynomiality and positivity. 
Beyond leading logarithmic order one would need to include the non-pointlike hadronic contri- 
bution. Note that the same process in a different kinematic region, namely at low energy and 
high squared momentum transfer —t ~ Q^, gives information on the generalized distribution 



a quark- ant iquark pair to the 



amplitudes (GDA) of the photon, which describes the coupling o: 
two photons and are connected to the photon GPDs by crossing [5]. DVCS process on a proton 
target has been analyzed in detail theoretically and are also being accessed in experiments jo]. 
The proton GPDs are richer in content than the ordinary parton distributions (pdfs). In the 
forward limit of zero momentum transfer they reduce to pdfs and their x moments give nucleon 



2 



form factors. An interesting physical interpretation of GPDs has been obtained in [3] by taking 
their Fourier transform with respect to the transverse momentum transfer. When the longitu- 
dinal momentum transfer is zero, this gives the distribution of partons in the nucleon in the 
transverse plane. They are called impact parameter dependent parton distributions (ipdpdfs) 
q{x, 6"*"). In fact they obey certain positivity constraints which justify their physical interpreta- 
tion as probability densities. This interpretation holds in the infinite momentum frame (even 
the forward pdfs have a probabilistic interpretation only in this frame) and there is no relativis- 
tic correction to this identification because in light-front formalism, as well as in the infinite 
momentum frame, the transverse boosts act like non-relativistic Galilean boosts. When the 
nucleon is transversely polarized, the unpolarized impact parameter dependent pdf is distorted 
in the transverse plane. A combination of chiral odd GPDs in impact parameter space gives 
information on the correlation between the spin and orbital angular momentum of the quarks 
inside the target jsf. Fourier transform (FT) with respect to the skewness ( gives rise to a 
diffraction-like pattern jo]. Thus the GPDs in effect give a complete (Lorentz invariant) 3 D 
picture of the proton in position space. While the proton is known to be a composite particle, 
it is interesting to access the partonic structure of the photon probed in high energy processes. 
As the proton GPDs are richer in content than the ordinary pdfs, photon GPDs can shed more 
light on the partonic content of the photon. 

In this letter, we calculate the photon GPDs using overlaps of light-front wave functions. 
We take the momentum transfer to be purely in the transverse direction, unlike where 
the momentum transfer was taken purely in the lightcone (plus) direction. We keep leading 
logarithmic terms and the mass terms coming from the vertex; and upto leading order in 
electromagnetic coupling and zeroth order in strong coupling. There are contributions only 
from the diagonal (particle number conserving) overlaps. When there is nonzero momentum 
transfer in the longitudinal direction, there are off-diagonal particle number changing overlaps 



as well, similar to the proton GPDs 



10| . The two particle light-front wave functions of the 



photon can be calculated analytically using perturbation theory. Taking a Fourier transform 
with respect to the momentum transfered in the transverse direction. A"*-, we express the GPDs 
in the transverse impact parameter space. 
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GPDs of the photon 



The GPDs for the photon can be expressed as the following off-forward matrix elements 
defined for the real photon (target) state {4]: 



= 



^e^(7(P') I miMy-) I 7(P)); 

Stt 

^e^(7(P') I V^(0)7V^(l/-) I liP)). 

OTT 



(1) 



F'^ contributes when the photon is unpolarized and F'^ is the contribution from the polarized 
photon. We have chosen the light-front gauge = 0. As pointed out in [4], there is also 
the photon operator which mixes with the quark operator, it contributes at the same order in 
a in the scattering amplitude. However here we calculate the matrix element rather than the 
amplitude of the process and the contribution to the matrix element of the photon operator 
comes at zeroth order in a. P^ and P'' can be calculated using the Fock space expansion of the 
state, which can be written as 



7( 



0'1,(J2 



(P)) = ViV a^{P, X)\0)+J2 / {^^1} / {dk2} ^2{27r fP+6^{P - h - h 



(j)2{ki,k2,ai,a2)b\ki,ai)d\k2,cr2) \ 0) 



(2) 



where vA^ is the overall normalization of the state; which in our calculation we can take as 



unity as any correction to it contributes at higher order in a. {dk} = J 



dk+d?k^ 



^2 is the 



two-particle {qq) light-front wave function (LFWF) and ai and (72 are the helicities of the quark 
and antiquark. The wave function can be expressed in terms of Jacobi momenta Xi = ^ and 
= k^ — XiP-^. These obey the relations ^jXj = 1, = 0- The boost invariant LFWFs 

are given by il)2{xi.,ql) = 4>2^ ■ '^2{xi-iqi) can be calculated order by order in perturbation 
theory. The two-particle LFWFs are given by |l2| 

1 



^2 _ rn?+{>l^) 
x{l—x) 

\ — X 



2(27r) 



[a^ ■ q' 



-a 



X 



m 



x(l 



-a 



X] 



(3) 
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where we have used the two-component formahsm ll|, ll2| and m is the mass of q{q). A is 
the hehcity of the photon and si, S2 are the hehcities of the q and q respectively. There is no 
contribution to the matrix elements that we consider from the single particle sector of the Fock 
space expansion. The leading term is the two-particle contribution, which can be written as, 



(fq^dxiS{x — Xi)^p2i^iy ^Z"*" "~ (1 "~ Xi)A^)ip2{xi, q^) 



(4) 



Here we have suppressed the helicity indices and the sum over them. The momentum transfered 
square t is given by t = {P — P'Y = ~{^^)^ ■ The first term is the contribution from the 
quarks and the second is the contribution from the antiquark in the photon. As the light-cone 
momentum fraction Xi has to be always greater than zero, the first term contributes when 
1 > X > and the second term for — 1 < x < 0. Using the LFWFs each component can be 
calculated separately. We calculate in the same reference frame as j9l] . Note that the light cone 
plus momentum of the target photon is non-zero. Finally we get for the unpolarized photon 

2 

= [((1 - xf + x'){h + I2 + Lh) + 2m'h]e{x)e{l - x) 

cac^ r 

[((1 + + + ^2 + L'l!,) + 2m'l',\ e{-x)e{l + x) (5) 



Here the sum indicates sum over different quark flavors; L = —2m'^+2m'^x{l—x) — {A^) (1— x) 
L' = —2m?' — 2m?x{l + x) — (A^)^(l + x)^; the integrals can be written as, 

^ r A' 



h 



D 



nLog 



.fi"^ — m2x(l — x) + m^. 
da 



DD' 



IT 







P(x,a, (A- 



(6) 



where D = (g^)^ — m^x(l — x) and D' = (g-*-)^ + (A^)^(l 
x) + m^, and P(x, a, (A-*-)^) 



-X 



2g^- A^(l-x)-m2x(l 



— m^x(l — x) + + a(l — a)(l — x)^(A-'-) . At zeroth order in 
as the results are scale dependent, this scale dependence in our approach comes from the upper 
limit of the transverse momentum integration A = Q. /i is a lower cutoff on the transverse 
momentum, which can be taken to zero as long as the quark mass is nonzero. Leading order 
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evolution of the photon GPDs has been calculated in for non-zero (. The mass terms in the 
vertex eive subdominant contributions which we included. 





FIG. 1: (Color online) (a) Plot of unpolarized GPD F'^ vs x for fixed values of —t in GeV'^ and (b) 
polarized GPD Fi vs x for fixed values of —t in GeV'^^ A = 20GeV. In both plots the normalization 
factor is chosen to compare with when t = 0. 



For the antiquark contributions we have similar integrals 

cPq"^ . r A' 



H 



TiLog 



HH' 



da- 



77 



(7) 



where H = {q^f + m'^x{l + x) + and H' = {q^f + (A-L)^(I + xf + 2q^ ■ + x) + 

m?x{l + x) + m^, and Q{x, a, (A-^)^) = m^x(l + x) + + a(l — a)(l + x)^(A^)^. 
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FIG. 2: (Color online) (a)Plot of impact parameter dependent pdf q{x,b) vs x for fixed b values and 
(b) q{x, h) vs b for fixed values of x where we have taken A = 20 GeV and lS.max= 3 GeV where A^ax 
is the upper limit in the A integration, b is in GeV^^ and b) is in GeV^ . 



For polarized photon the GPD F'^ can be calculated from the terms of the form e|ej^* — e^e|* 
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J]. We consider the terms where the photon hehcity is not flipped. This can be written as, 



= 5^^[(a^'-(l-^)')(/i + /2 + i:/3) + 2m%]^(x)0(l-x) 



ryp, 

^ ^ (x^ - (1 + + + L'l',) + 2mH', e{-x)e{l + X 



(8) 
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FIG. 3: (Color online) (a)Plot of impact parameter dependent pdfs q{x, b) vs x for fixed b values and 
(b) q{x, b) vs b for fixed values of x where we have taken A = 20 GeV and Amax= 3 GeV where Amax 
is the upper limit in the A integration, b is in GeV^^ and q{x, b) is in GeV^. 
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In analogy with the impact parameter dependent parton distribution of the proton, we 
introduce the same for the photon. By taking a Fourier transform with respect to the transverse 
momentum transfer A"*- we get the GPDs in the transverse impact parameter space. 



q{x, 6"* 



1 



(2n) 



^ / ArfAJo(A6)F« 
2n J 



(9) 



q{x, 6"* 



(27r)2 
1 



27r 

: lA^I and b 



AdAJo{Ab)F'J; 



(10) 



16 1. In the numerical calculation, we 



where Jo{z) is the Bessel function; A 
have introduced a maximum limit A^ax of the A integration which we restrict to satisfy the 
kinematics — t << js, 9, Is], Q- q{x,b-^) gives the distribution of partons in this case 
inside the photon in the transverse plane. Like the proton, this interpretation holds in the 
infinite momentum frame and there is no relativistic correction to this identification because in 
light-front formalism, as well as in the infinite momentum frame, the transverse boosts act like 
non-relativistic Galilean boosts. q{x, b^) gives simultaneous information about the longitudinal 
momentum fraction x and the transverse distance b of the parton from the center of the photon 
and thus gives a new insight to the internal structure of the photon. The impact parameter 
distribution for a polarized photon is given by q{x, b^). 



Numerical results 

We have plotted the unpolarized GPD and the polarized GPD F"^ for the photon in 
Figs 1(a) and (b) respectively as functions of x and for different values of t = — (A-*-)^. In 
all plots we took the momentum transfer to be purely in the transverse direction. We took 
m = 3.3 MeV; A = Q = 20GeV, /i = MeV and Amax= 3.0 GeV where Amax is the upper 
limit of the A-*- integration in the Fourier transform. It is to be noted that the photon structure 
function F^(x,(5^) has been explored over a wide kinematical range, namely 0.001 < x < 0.9 
and 1.9 < < 780GeV^ However, here we restrict ourselves to study the general features 
of the photon GPDs at a fixed scale rather than the scale evolution. We have divided the GPDs 
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FIG. 4: (Color online) Plots of impact parameter dependent pdfs (a) q{x,b) and (b) q{x,b) vs b for 
a fixed x; where we have taken A = 20 GeV for different values of A^ax in GeV. b is in GeV~^ and 
q{x,b) as well as q{x,b) are in GeV^. 

by the normaUzation constant to compare with j4| in the hmit of zero t. Indeed they agree. As 
—t increases, F'^ becomes more and more asymmetric with respect to a; = 1/2 : this asymmetry 
is prominent for lower values of x which is expected as the A-*- or t dependence is associated 
with a (1 — x)^ factor (see the analytic expressions). The slope of the polarized GPD changes 
with increasing (A-*-)^. Note that in {4] as well as in the solid line in figs 1(a) and 1(b), the 
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subleading mass terms are not taken into account. It is to be noted that in all plots we have 
taken < x < 1 for which the contribution comes from the active quark in the photon {qq). 
As X — )■ 1, most of the momentum is carried by the quark in the photon and the GPDs become 
independent of t. The Fourier transform (FT) of the unpolarized GPD is plotted in Fig. 2. 
Fig. 2 (a) shows the plot of the impact parameter dependent pdf of the photon as a function 
of X and for fixed impact parameter b =\ b-^ \. Fig. 2 (b) shows the same but as a function 
of b and for fixed x. The smearing in b-^ space reveals the partonic substructure of the photon 
and its 'shape' in transverse space. In the ideal definition the Fourier transform over A should 
be from to oo. In this case the independent terms in F'^ and F'^ would give 5^(6^) in 
the impact parameter space. This means in the case of no transverse momentum transfer, the 
photon behaves like a point particle in transverse position space. The distribution in transverse 
space is a unique feature accessible only when there is non-zero momentum transfer in the 
transverse direction. From the plots it can be seen that for fixed b, q{x, b) decreases slowly 
with X till X ~ 1/2 then increases. The behavior in impact parameter space is qualitatively 
different than a dressed quark target and also from phenomenological models of proton GPDs. 
For a dressed quark target the leading contribution to the GPD comes from the single particle 
sector of the Fock space expansion, which in impact parameter space gives a delta function 
peak. This contributes at x = 1. For large values of x the peak in impact parameter space 
becomes sharper and narrower, that means there is a higher probability of finding the active 
quark near the transverse center of momentum [l^. In the phenomenological parametrization 
of proton GPDs where a spectator model with Regge-type modification was used, the GPDs 
have a different behavior in impact parameter space, the u quark GPDs increase with increasing 
X for fixed 6, reaches a maximum, then decrease. The peak decreases with increasing b jl^ . 
In the case of a photon there is no single particle contribution, and the distribution in b space 
purely reveals the internal qq structure of the photon. Here near x ~ 1/2 the peak in b space is 
very broad which means that the parton distribution is more dispersed when the q and q share 
almost equal momenta. The parton distribution is sharper both for smaller x and larger x. In 
Figs 3 (a) and (b) we have plotted the polarized distribution in impact parameter space. Fig. 3 
(a) shows it as a function of x for fixed b values and fig 3 (b) shows it as a function of b for fixed 
X values. The slope decreases for higher b. The sign of the GPD changes at x = 1/2, at which 
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point the GPD and the pdf in impact parameter space becomes zero. The distributions are 
approximately symmetric about a; = 1/2 in impact parameter space; as x — 1 the distribution 
increases sharply as the GPDs become independent of t; this is similar to q{x,b). For fixed x, 
q{x, b) as a function of b becomes broader as x increases until x = 1/2. For larger values of x, 
it changes sign. In the plots we have taken the upper limit of the Fourier transform to be much 
smaller than A. The dependence of q{x, b) and q{x, b) on Amax is shown in Figs. 4 (a) and (b) 
respectively. As Amax increases, the distribution becomes sharper in impact parameter space. 
This shows that larger momentum transfer probes the partons near the transverse center of the 
photon. 

Conclusion 

We presented a first calculation of the generalized parton distributions of the photon, both 
polarized and unpolarized, when the momentum transfer in the transverse direction is non- 
zero; at zeroth order in and leading order in a; we calculated at leading logarithmic order 
and also kept the mass terms at the vertex. We took the skewness to be zero. We express 
the GPDs in terms of overlaps of the photon light-front wave functions. We considered the 
matrix elements when the photon helicity is not fiipped. When the momentum transfer in the 
transverse direction is non-zero, one also has helcity fiip contributions, which will be treated 
in a later work. In our case only the diagonal parton number conserving overlaps contribute. 
We considered both the quark and the antiquark contributions. The GPDs thus probe the two 
particle qq structure of the photon. Taking a Fourier transform (FT) with respect to A_l we 
obtain impact parameter dependent parton distribution of the photon. We plot them for both 
polarized and unpolarized photon. The parton distributions in impact parameter space show 
distinctive features compared to the proton and also compared to a dressed quark, which can 
be taken as an example of a spin 1/2 composite relativistic system consisting of a quark and a 
gluon. It is to be noted that a complete understanding of the photon GPDs beyond leading logs 
would require also the non-pointlike hadronic contributions which will be model dependent ^. 
However, the GPDs of the photon calculated here may act as interesting tools to understand 
the partonic substructure of the photon. Accessing them in experiment is a challenge. On 
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the theoretical side, the next step would be to investigate the photon GPDs when there is 
non-zero momentum transfer both in the transverse and in the longitudinal direction as well 
as a perturbative study of the general properties of GPDs like positivity and polynomiality 
conditions and sum rules. 
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